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Although materials that exhibit nearest-neighbor-only antiferromagnetic interactions and geomet- 
rical frustration theoretically should not magnetically order in the absence of disorder, few such sys- 
tems have been observed experimentally. One such system appears to be the pyrochlore Tb2Ti207. 
However, previous structural studies indicated that Tb2Ti207 is an imperfect pyrochlore. To clarify 
the situation, we performed neutron powder diffraction (NPD) and x-ray absorption fine structure 
(XAFS) measurements on samples that were prepared identically to those that show no magnetic 
order. The NPD measurements show that the long-range structure of Tb2Ti2 07 is well ordered 
with no structural transitions between 4.5 and 600 K. In particular, mean-squared displacements 
(tt 2 's) for each site follow a Debye model with no offsets. No evidence for Tb/Ti site interchange 
was observed within an upper limit of 2%. Likewise, no excess or deficiency in the oxygen stoi- 
chiometry was observed, within an upper limit of 2% of the nominal pyrochlore value. Tb Lm- 
and Ti K-edge XAFS measurements from 20-300 K similarly indicate a well-ordered local structure. 
Other aspects of the structure are considered. We conclude that Tb2Ti207 has, within experimental 
error, an ideal, disorder-free pyrochlore lattice, thereby allowing the system to remain in a dynamic, 
frustrated spin state to the lowest observed temperatures. 

PACS numbers: 75.25.+Z, 75.40.Gb, 75.50.Ee, 75.50.Lk, 61.12.Ld, 61.10.Ht 



I. INTRODUCTION 

The term geometrical frustration*'^ indicates when 
the spins on a well-ordered lattice interact magnetically, 
yet individual interactions cannot reach their minimum 
energy configuration because of competing magnetic in- 
teractions from other sites. A simple example occurs on 
a two-dimensional, corner-shared triangular lattice with 
antiferromagnetic nearest-neighbor coupling, otherwise 
known as a kagome lattice. The three dimensional ana- 
logue to such a lattice is a corner-shared tetrahcdral lat- 
tice, and is realized in the pyrochlore systems, A2-B2O7. 

Theoretically, it has been shown that Heisenberg, 
antiferromagnetic nearest-neighbor interacting spin sys- 
tems^ do not possess a transition into an ordered state. 
This makes corner-sharing tetrahedral spin systems with 
antiferromagnetic exchange ideal candidates for three- 
dimensional, low temperature spin-liquids. A spin-liquid 
state may not form if further interactions such as crys- 
tal field anisotropy, 5 , long-range exchange or dipolar in- 
teractions perturb the classical spin models^ However, 
since any amount of lattice disorder can precipitate a 
spin-glass phase, 7 geometrically-frustrated systems that 
do not magnetically order to the lowest temperatures are 
extremely rare£ This situation has inhibited experimen- 
tal studies that attempt to isolate the effects of frustra- 
tion. For instance, lattice imperfections in several geo- 
metrically frustrated systems are thought to relieve the 
frustration allowing the system to freeze. Neel order 
has been observed in the Jarosites, KCr3(OD)e(S04)2, 
which contains Cr 3+ , S=| ions on a kagome lattice^ 
due to vacancies on the magnetic site. Site disorder 
in SrCr9 P Gai2-9pOi9fi& for example, facilitates the for- 



mation of a spin glass state, as does intersite mixing in 
Gd3Ga 5 12 ii^ 

Significantly, the pyrochlore lattice^ appears to 
be capable of producing a spin-liquid state. The 
magnetic-oxide pyrochlores have the chemical formula 
A2-E>20(1)60(2) (space group Fd3m). Here, we focus 
on the case where the A site is occupied by a tri- 
valent rare-earth ion with eightfold oxygen coordination, 
ylO(2) 2 0(l)6, and the B site is occupied by a tetra-valent 
transition metal ion with sixfold oxygen co-ordination, 
50(1)6- The A and B sites individually form infinite 
interpenetrating sublattices of corner-sharing tetrahedra. 
A high degree of magnetic frustration is exhibited on this 
lattice when either site is occupied by antiferromagneti- 
cally coupled ions. 

Materials with the A2-B2O7 formula unit can have var- 
ious structures, some of which can be thought of as de- 
fected versions of the fluorite (CaF2) unit cell. The cu- 
bic fluorite cell has Ca atoms at the center of eight F 
atoms at the corners of a cube. This cell is analogous 
to one eighth that of the oxide pyrochlore with a ran- 
dom distribution of the metal cations, one oxygen va- 
cancy (the 8a site), and the positional parameter x for 
the oxygens on the 48/ site equal to 3/8. The rare-earth 
titanate pyrochlores have x ~ 0.328. 1 — Obviously, these 
structures are related, but the mixing of the metal ions 
in the fluorite structure creates multiple exchange con- 
stants, usually resulting in a static magnetic structure. 
Other materials with the stoichiometry A^B^Oi include 
the weberites, which are closely related to the pyrochlore 
structure, but only half the i?-site ions have the sixfold 
oxygen coordination, while the others are linked through 
only four vertices. 
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Despite the inherent geometrical frustration in the 
ideal pyrochlore lattice, most pyrochlores order magnet- 
ically. Tb2Ti207, on the other hand, appears to be the 
rare case of a pyrochlore that does not order magnet- 
ically, exhibiting a high degree of frustration (the large 
Tb 3+ moment has a Curie- Weiss temperature of ~-19 K) 
down to 70 mKpiSiiSiil However, the pyrochlore system is 
also capable of producing a lattice with fairly subtle dis- 
order that can precipitate a spin-glass phase, such as in 
the Mo-Mo near-neighbor pairs in Y2Mo2C , 7ii2i 1 2 A sim- 
ilar mechanism may also precipitate the spin-glass phase 
in Tb2Mo207i2£ In fact, the lack of such frozen magnetic 
states in Tb2Ti20y is somewhat surprising, given the re- 
sults of previous structural measurements. The crystal 
structure of T^T^Oy was first investigated by Brixnei«2i 
in 1964. The most detailed structural study is reported 
by van de Velde et alw^ The International Centre for 
Diffraction Data quotes van de Velde's data as the stan- 
dard for Tb2Ti207. Their x-ray diffraction study finds 
that the measured Tb2Ti2C>7 sample has approximately 
10% disorder on the metal cation and oxygen sublattices, 
and the authors conclude that Tb2Ti20y is an imperfect 
pyrochlore. The study we report below was motivated by 
this result and attempts to definitively answer the ques- 
tion of lattice imperfections in Tb2Ti2C>7, a significant 
parameter in any model that would successfully describe 
the magnetic nature of Tb2Ti2 07. 

The present work reports neutron powder diffraction 
(NPD) measurements of the average structure, and care- 
fully considers the thermal dependence of the mean- 
squared displacement parameters as a means of sepa- 
rating the non-Debye (presumably static) component of 
any positional disorder. In addition, the x-ray absorption 
fine-structure (XAFS) technique is used to check for dis- 
order in the local environment around Tb and Ti atoms, 
and in particular for the kind of -B-site correlated disor- 
der that occurs in Y2M02O7. Finally, the results from 
the two techniques are combined to determine to what 
degree correlations exist between several atom pairs. 

The rest of this paper is organized as follows: Sample 
preparation and experimental details are in Sec. ^ Data 
analysis procedures and results are reported in Sec. IIIII 
Consequences of these results, comparisons to previous 
measurements, and further discussion are in Sec. II VI 
and the final conclusions are in Sec. [V] 
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FIG. 1: Neutron diffraction data collected at IPNS from pow- 
der Tb2Ti2 07 at 45 K as a function of the momentum transfer 
q. The solid curve is the best fit from the Rietveld refinement 
with the Fd3m structure. The bottom curve shows the dif- 
ference between the measured and calculated intensities. 



NRU reactor at Chalk River Laboratories. The ground 
powder was packed in a vanadium can and cooled in 
a displex (IPNS) or a helium-bath cryostat (NRU). A 
Si(113) crystal was used at the NRU reactor to produce 
a monochromatic neutron beam with a wavelength of 
1.3283 A. The large difference in the scattering lengths, 
in particular the negative scattering length of titanium, 
creates a large contrast with the terbium scattering. This 
contrast allows for an accurate determination of the crys- 
tal structure properties, and especially for any possible 
cation antisite disorder. 

For XAFS measurements, the polycrystalline 
Tb2Ti2C>7 sample was re-ground and passed through 
a 20 /im sieve. The sieved powder was uniformly dis- 
tributed over adhesive tape, cut into strips and stacked 
to obtain absorption edge steps of 0.8 absorption lengths 
at the Ti K edge and 1.0 absorption lengths at the Tb 
Liu edge. Samples were then placed into a LHe-flow 
cryostat. Transmission XAFS data were collected at the 
Stanford Synchrotron Radiation Laboratory (SSRL). Tb 
Lin-edge (7515 eV) data were collected on beamline 4-3 
with a 1/2-tuned Si(220) double-crystal monochromator. 
Ti K-edge (4965 eV) data were collected on beamline 2-3 
with a 1/2-tuned Si(lll) double-crystal monochromator. 



II. EXPERIMENTAL DETAILS 



III. RESULTS 



Polycrystalline Tb2Ti207 was prepared by a conven- 
tional solid state reaction in the same manner as in Rcf. 
[l5l namely, by firing stoichiometric amounts of Tb407 
and TiC>2 at 1350° for several days with intermittent 
grindings to ensure a complete reaction. The neutron 
diffraction measurements were performed on the time-of- 
flight Special Environment Powder Diffractometer at the 
Intense Pulsed Neutron Source (IPNS) at Argonne Na- 
tional Laboratory and on the C2 diffractometer at the 



The rare earth oxide pyrochlores crystallize in the 
space group Fd3m as mentioned above, with 4 crystal- 
lographically independent sites (A 3+ in position I6d at 
(|, |, |), B 4 + in 16c at (0, 0, 0), 0(1) in 48/ at (x, §, §) 
and 0(2) in 8b at (§, §, §)±£). The position of the 48/ 
site is defined by only one variable, x, which is sensitive 
to lattice effects, including the absolute size and the rel- 
ative difference in size of the metal ions, as well as to 
various forms of disorder^ 
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FIG. 2: Mean-squared displacement parameters, 
4 crystallographic sites of Tb2Ti2 07 as a function of temper- 
ature are shown in panels (a)-(c). The data were collected at 
the IPNS (closed symbol) and the NRU reactor (open sym- 
bol). Solid lines are best fits using a Debye model (see Table 
EJ. Panels (d) and (e) show the temperature dependence of 
the positional parameter x of the O(l) oxygen (48/ site) and 
the lattice parameter a, respectively. 



The neutron diffraction data were analyzed with a Ri- 
etveld method using the General Structure Analysis Sys- 
tem (GSAS) refinement packaged Figure Q shows the 
45 K data together with the best refinement and the 
residual, shown at the bottom. No additional peaks, as- 
sociated with either a second phase or a magnetic unit 
cell, were observed over the whole temperature range 
from 4.5-600 K. However as previously reported^ an 
anomalous background due to short range magnetic cor- 
relations was observed at the lowest temperatures. The 
atomic displacement parameters of all atoms were as- 
sumed to be isotropic during the fitting process since no 
improvement in the refinement was obtained by using 
anisotropic parameters. This is in contrast to Y2M02O7, 
where anisotropic parameters for the A site were deemed 
necessary^ The results of the refinement are summa- 
rized in Tabic HJ The possibility of antisite disorder (site 
interchange between Tb and Ti atoms) was also consid- 
ered, such as could occur if the system had a propensity 
toward a fluorite-like structure i^SiSi We find no evidence 
for any such disorder, and place an upper limit of 2% 
on metal-site mixing. Similarly, we performed fits al- 



TABLE I: NPD fit results to IPNS data from Tb 2 Ti 2 7 us- 
ing the Fd3m space group. Reported values are from fits at 
45 K, except for the Debye fit results, which are over the full 
temperature range of the data (45 - 600 K). Static disorders 
were obtained with the Debye model fit shown in Fig. [5] All 
site occupancies are fixed at unity. 



General fit characteristics: 
Banks included 
Total data points 
Total measured reflections 
# of variables 

reduced \ 2 

Rp(%) 

wRp(%) 



± 145°, ±90° 

8792 

362 

9+4 for background 



1.802 

6.90 

4.51 



Atom x 



V 



,(A) 6 D (K) uLtic(A 2 ) 



Tb 
Ti 
0(1) 
0(2) 



1/2 


0.3285(1) 
3/8 



1/2 


1/8 
3/8 



1/2 


1/8 
3/8 



0.0018(2) 
0.0028(4) 
0.0039(2) 
0.0033(3) 



220(10) 
610(50) 
760(30) 
980(80) 



0.0006(7) 
0.0015(12) 
0.0010(12) 
0.0008(12) 



a 10.13315(12) A 



Tb-0(2)-Tb 109.47(1)° 

Tb-0(1)-Tb 106.41(3)° 

Tb-0(1)-Ti 106.41(2)° 

Ti-0(1)-Ti 132.12(5)° 



lowing the oxygen occupancy fractions to vary, includ- 
ing the possibility of oxygen on the nominally vacant 8a 
site. These measurements provide no evidence of excess 
or deficient oxygen compared to the nominal pyrochlorc 
structure within an upper limit of 2%. Both of these lim- 
its are considerably lower than the disorder reported by 
van de Velde et al^ 

Figure |21 shows the mean-squared displacements, u 2 , 
for all four crystallographic sites and the positional pa- 
rameter, x, as a function of temperature. No structural 
anomalies, within uncertainty, were observed. The it 2 's 
for each site were fit with a thermal model, including a 
temperature-independent offset: 



U meas(-^) — ""static + ^thermal (O * 



(1) 



The temperature-dependent part of the mean-squared 
displacement parameters, u 2 hermal (T) , is given in general 
by a collection of quantum harmonic oscillators: 



^thermal CO n 



p(co)coth 



2k B T 



- (2) 

U! 



where m is the atomic mass and p is the phonon den- 
sity of states at the given site. In the Debye model, the 
density of states is empty above a cutoff frequency, wd, 
and quadratic below it: p — Zuj 2 /w D . The Debye tem- 
perature is then just 9d = Rwd/^b- Note that zero-point 
motion exists such that even in the absence of a static 
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FIG. 3: XAFS data at 20 K measured above the (a) Tb L 3 - 
edge and (b) Ti K-edge as a function of photoelectron wave 
number, k. 



offset, /thermal ls n °t zero a ^ T=0, and is, in fact: 



thermal 



(T = 0) 



3?r 



4?7lfcB0D 



The temperature-independent offset in addition to this 
zero-point motion, u 2 tatic , when observed, can be a strong 
indication of the presence of static (non-thermal) disorder 
in a system. The temperature dependence of the u 2 pa- 
rameters is described well by the Debye model as shown 
by the fits in Fig. El In Table [Q we summarize these fits. 
The itg tatic parameters are negligibly small at all atomic 
sites, and the Debye temperatures are appropriate for the 
elements they represent. 

Investigations of the local environment around the 
metal ions were performed using temperature depen- 
dent XAFS above the Tb Lm edge and the Ti K edge. 
XAFS oscillations occur above a core-level x-ray absorp- 
tion edge, and are due to the interaction between the 
outgoing photoelectron from the absorbing atom and the 
component of this photoelectron that is back-scattered 
by neighboring atoms. The oscillations are periodic in 
the pair distance, r, and the photoelectron wave num- 
ber, k = y / 2m(E — Eo)/h, where E is the incident x- 
ray energy and Eq is the energy at the absorption edge. 
XAFS analysis gives information about the radial dis- 
tance distribution around the absorbing atomic species. 
In particular, pair distances, r,-, and pair-distance dis- 
tribution widths, erf, can be obtained for the first few 
atomic shells, i, out to about 5 A. These XAFS data 
were analyzed using the UWXAFS package^. First, the 
XAFS function \ = M/Mo — 1 is extracted after deter- 
mining the atomic background absorption, /x , obtained 
with the UWXAFS program AUTOBK by passing a cu- 
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FIG. 4: Magnitude of the Fourier transformed (FT) XAFS 
data from Fig. El (thick-dotted lines) and the real part of 
the XAFS (thin-dotted lines) as function of distance from 
the probe atoms, (a) Tb L3-edge and (b) Ti K-edge. XAFS 
data in the range of 4.0 — 13.0 A -1 are used for the Fourier 
transform with a 0.5 A -1 wide Hanning window. Solid lines 
are the best fits and the vertical lines indicate the fit regions. 



bic spline through the absorption data, fi. The fc-space 
XAFS data for both measured edges are shown in Fig. 
Multiple scans were obtained at each edge and tempera- 
ture, with excellent reproducibility. 

Figure 01 shows the magnitude and the real part of the 
Fourier transformed XAFS for both edges. Note that 
the peaks are shifted on the f axis from their true bond 
lengths due to the phase shift of the back-scattered pho- 
toelectron. This phase shift, as well as other details of the 
complex back-scattering function, differ for each atomic 
species, although they are well approximated by the the- 
oretical treatment provided by the FEFF8 code.? 9 De- 
tailed fits are therefore necessary to obtain quantitative 
information. Solid lines show the best fit of a multiple 
shell model that assumes full occupancy of each of the 
sites, but allows the and the a 2 parameters for each 
shell to vary, with some exceptions. In particular, the 
Tb-Ti and Tb-Tb pairs near 3.60 A are constrained to 
have the same pair distance, as are the Ti-Ti and Ti- 
Tb pair distances. The fit results are given in Table ITT1 
and Fig. Additional details of the XAFS analysis and 
the determination of the uncertainty of parameters have 
been given elsewhere iSS^li 

The a 2 data were analyzed analogously to the NPD 
u 2 parameters described above, except that a correlated- 
Debye model^S^ was used to account for correlated mo- 
tion/displacements of the atom pairs. This model for a 2 
is the same as the Debye model (Eq. ^and|2) except that 
a correlated-Debye temperature, Q c d, is used to differ- 
entiate from the regular Debye model, the reduced mass, 
(i, of the atom pairs is used, and the phonon density of 
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TABLE II: XAFS fit results at 20 K. The overall amplitude reduction factors, Sq, were determined to be 0.96(9) for the Tb 
Lin-edge fits and 0.80(6) for the Ti _ft"-edge fits. The number of fit parameters (10 and 13, respectively), is much less than 
the number of independent data points determined by Stern's rule— (19 for each). TV is the coordination number for the fully 
occupied model. Bond lengths of atomic pairs determined using XAFS (r) are compared to the NPD results (t*npd) at 45 K. 
The correlated-Debye temperatures, c d, and the temperature- independent offsets, u% atil ., were determined by fitting with a 



correlated-Debye 


mode&S 


over the measured temperature range (20-300 K). The correlation parameter cf. 


is given by Eq. 2] 


Atomic pair 


N 


r(A) 


»"npd(A) 


e cD (K) 


^statici-^ ) 




Tb-0(2) 
Tb-O(l) 
Tb-Ti 1 
Tb-Tb 1 


2 
6 
6 
6 


2.199(6) 
2.490(9) 
3.598(6) 
3.598(6) 


2.1939(6) 
2.4958(7) 
3.583(1) 
3.583(1) 


875(170) 
450(20) 
300(20) 
269(11) 


-0.0002(7) 
0.0004(5) 
-0.0004(8) 
-0.0003(3) 


0.61(7) 
0.18(7) 
0.40(15) 
0.65(9) 


Ti-O(l) 
Ti-Ti 2 
Ti-Tb 2 


6 
6 
6 


1.985(7) 
3.596(7) 
3.596(7) 


1.9699(5) 
3.583(1) 
3.583(1) 


600(50) 
500(60) 
300(15) 


-0.0009(10) 
-0.0005(7) 
-0.0003(5) 


0.60(7) 
0.55(8) 
0.38(10) 



1 Tb— Ti & Tb— Tb bond length constrained together 

2 Ti— Ti & Ti— Tb bond length constrained together 




FIG. 5: Debye- Waller factor(<r 2 ) as a function of temperature, 
(a) Tb-Ti and Ti-Tb pairs, and (b) Tb-O(l), Tb-0(2) 
and Ti— O(l) pairs determined at Tb L^-edge and Ti K-edge. 
Solid lines are the best fits with the correlated-Debye model 
(see Table ED. Ti-Ti pairs are omitted for clarity. 



states for a given pair at distance R is given by 



P 



,3 



D 



1 



sin{ijjR/c) 



ujR/c 



(3) 



where lo-q is the usual Debye cut-off frequency and 
c = uji)/k-£,. The Debye wave vector is given by &d = 
(6tt 2 /V) 1/3 , where V is the mean volume per atom in 
the material. The expression in brackets of Eq. Stakes 
into account the correlated motion of the atom pair. This 
model fits the temperature dependence of the measured 
<7 2 's well (Fig. [SJ with only negligible static displace- 



ments, cr 2 tatic 's, as shown in Table ITT1 This lack of pair- 
distance disorder agrees well with the lack of positional 
disorder determined from the NPD results above. 

A Tb/Ti site interchange model was also considered 
for these XAFS data. However, XAFS is insensitive to 
such site interchange in this system because correlations 
in the fits between the (nominally same) Tb-Tb and Tb- 
Ti pair distances (as well as the Ti-Tb, Ti-Ti distances), 
their cr 2 's, and the alleged interchange render the results 
inconclusive. 



IV. DISCUSSION 

The primary result of the above measurements is that 
the Tb2Ti207 lattice is extremely well ordered, both 
on average (NPD) and locally (XAFS). This result is 
supported by the high quality of the fits, the low es- 
timated errors, and the temperature dependence of the 
mean-squared displacement parameters. These data thus 
strongly support the conclusion made by previous au- 
thors that Tb2Ti2C>7 is an excellent laboratory to study 
the effects of geometric frustration in the absence of lat- 
tice disorder. 

The discrepancy between our results and those of van 
de Velde et alM> must be attributable to either the 
different techniques used to measure the structure or 
to produce the sample. Since powder x-ray diffraction 
should be able to distinguish the cations, we must as- 
sume that the synthesis method used in van de Velde et 
al., namely, preparing Tb2Ti20y by the citrate method 
from TiO(OH)2 and a solution of Tb407 in nitric acid, re- 
sults in an imperfect pyrochlore lattice. This situation is 
not inconceivable, since considerably lower temperatures 
and shorter times are used in this method compared to 
the solid state reaction method used here^ 

Very recently, Lian et nlml reported an x-ray diffrac- 
tion study of the entire RE2Ti20y (RE=rare-earth) se- 
ries. Although they do not focus on Tb2Ti2 07 or on 
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the detailed search for lattice disorder presented here, 
they report a room-temperature x parameter for the 
oxygen 48/ site of 0.3281(5), which is consistent with 
our measurement of 0.3279(1), but not with the value 
of 0.317(2)24 found by van de Velde et alM A value of 
x w 0.328 is, however, consistent with that found in the 
RE2Ti207 series^ while a value less than 0.32 is a gross 
deviation. As noted earlier, this parameter is very sensi- 
tive to lattice distortions, which could be due to antisite 
disorder, oxygen excess or deficiency, etc. In fact, Min- 
ervini et alM find that the theoretical ideal value for x 
given no lattice disorder is ~0.3278, in excellent agree- 
ment with our measurements and those of Lian et al. 
This agreement, together with the fact that Lian et al. 
used a flux-growth synthesis technique and obtained sim- 
ilar results to those presented here, is therefore consistent 
with our interpretation that the main reason van de Velde 
et al. found Tb2Ti2 07 to be an imperfect pyrochlore is 
the synthesis method, not the measurement technique, 
nor is the disorder generic to the material itself. 

There are a few interesting points to note about the 
measured structure and its vibrational properties. One 
is that the single-site Debye temperatures (from NPD) 
do not scale as expected in a Debye solid, that is, by the 
square root of the ratio of their masses. In this way, one 
sees that the Tb atomic vibrations are governed by much 
softer phonons that the Ti atoms [a/ toti /wTb©D (Ti) « 
330 K, while 6 D (Tb)=220 K]. 

Another aspect to consider is the correlated displace- 
ments of the atom pairs. Since we have both neu- 
tron diffraction and XAFS measurements, we can ex- 
tract the correlations by combining the results from 
both techniques. The instantaneous distance (Atab) 
between two atoms (A, B) is equal to the difference of 
the instantaneous displacements of the individual sites, 
Ar A - Ar B . Taking a\ B =< Ar A 2 B >, u\ =< ArJ > 
and uaub^ = < ArU • Are >, we then have 

a AB = U A + U B ~ 2(j)U A UB. (4) 

The correlation coefficient, <j), describes the motion or 
displacements of an atom relative to the other atom, tak- 
ing the values +1 for two atoms moving together in the 
same direction, -1 for the atoms moving in fully opposite 
directions, and for uncorrelated motion/displacements. 
The </>'s of Tb2Ti2C>7 at the lowest temperature extracted 
from the XAFS and diffraction measurements are sum- 
marized in Tablelrfl The motions of the two 0(2) oxygens 
located at ~2.2 A from the Tb atom and the six 0(1) oxy- 
gens located at ~2.0 A from Ti atom are correlated in 
the motions to Tb and Ti, respectively, with <fi of ^60%. 
It is somewhat surprising that their motions are not more 
highly correlated with such short bond lengths, especially 
compared to, say, the copper-oxide superconductors. 35 
More interestingly, the Tb-O(l) pairs are nearly uncor- 
related in their motions, with a <fi of only 0.2. Therefore, 
the motions of the 0(1) atoms are governed mainly by 
the Ti sublattice, and not by the Tb sublattice. This 
result is consistent with the Tb sites exhibiting softer 



phonons than expected from the Ti vibrations. 

Given these vibrational properties, one should consider 
the charge distribution in the various ligands. One simple 
method is to determine whether the bond- valence sums^ 
for the Tb and Ti oxygen environments correspond to the 
expected values of +3.0 and +4.0, respectively. Indeed, 
we obtain +3.01 and +3.95 from the NPD data, indicat- 
ing that the charge is distributed as expected, given the 
measured coordination and bond lengths. Therefore, the 
measured lack of correlation between Tb and 0(1) atoms 
and the soft Tb phonons are most likely a generic feature 
of a well-ordered pyrochlore lattice, such as TI2MH2O7, 
or any of the other ferromagnetic pyrochlores. 

A natural question to ask then is "Why is Tb2Ti207 
structurally so well ordered?" To this end, we can only 
speculate as to why, for instance, Y2M02O7 exhibits dis- 
tortions while Tb2Ti207 does not. Ionic size arguments 
are marginally fruitful, since the ionic-size ratio^S for 
Tb2Ti207 (taItb = 1.72) is rather different than for 
Y2M02O7 [taItb = 1.57). However, both of these val- 
ues fall well within the range 1.46 to 1.80 typical of the 
pyrochloresii£ On the other hand, the negligibly small 
amount of static disorder on all atomic sites and be- 
tween all near-neighbor pairs might suggest that the in- 
teraction of corner sharing spins is not strong enough 
to distort the atomic sites, as has been conjectured for 
Y2Mo207ii2iiS This explanation is supported by the rel- 
atively small Curie- Weiss temperature of ~ —19 K in 
Tb 2 Ti 2 7 , compared to 200 K in Y 2 Mo 2 7 ^ 



V. CONCLUSION 

We have examined the bulk and local structure of 
the spin liquid pyrochlore Tb2Ti207 by powder neutron 
diffraction and XAFS at the Tb Lm and Ti K edges. 
These data indicate a high degree of crystalline order, 
suggesting that this compound forms a perfect pyrochlore 
lattice with no intersite mixing or anion disorder. Fur- 
thermore, analysis of the temperature dependence of the 
atomic displacement parameters and mean-squared dis- 
placements shows that the static lattice disorder is neg- 
ligible. Analysis of the Debye temperatures and the cor- 
relations between various site displacements show that, 
vibrationally, the Tb and Ti sublattices are partially de- 
coupled, and that the more plentiful 0(1) atoms are more 
tightly bound to the Ti sites. This compound, with 
its near absence of local and bulk structural disorder 
and antiferromagnetic nearest-neighbor interactionSfi^ 3 ^ 
is therefore an appropriate system to study the elusive 
3-dimensional, low-temperature, spin liquid. 
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